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PRE-SHUTTLE LIDAR SYSTEM RESK-\RCH 


INTRODUCTION 

This report describes the George Washington University ON CAMPUS 
technical support provided to NASA for the period April 1, 1984 
through July 31, 1986 under cooperative agreement NCCl-85. The effort 
was performed by Faculty and Students in the Electrical Engineering 
and Computer Science (EE&CS) Department coordinated with NASA Langlely 
Research Center personnel.* Included are the results of the initial 
phase of a simulation study in connection with Photomultiplier Tubes 
(PMTs) and associated networks and an analytical Investigation of the 
atmostpheric physics (Including multiscattering) leading to modeling 
studies in connection with Differential Absorption Lidar (DIAL) 
observations. This effort was in support of the ER-2 aircraft DIAL 
projects. Part I df this report describes the photomulitplier tubes 
simulation studies by Professor M. E. Zaghloul// and Part II of this 
report describes the analytical modeling studies of the atmospheric 
physics by Professor R. H. Lang//. 


*The authors are indebted to Robert Allen for his technical assistance 
and guidance. 

//Professors Lang and Zaghloul are members of the Electrical Engineering 
and Computer Science Department, George Washington University, 
Washington, D.C. 


PHOTOMULTIPLIER SIMULATION STUDIES 


by 


D.J. Rhee and M. E. Zaghloul 
The George Washington University 
Department of Electrical Engineering and Computer Science 
Washington, D.C. 20052 


ABSTRACT 

In this project a computer model to simulate the nonlienar 
response of the PMT due to the associated voltage-divider circuit is 
developed. An equivalent circuit which approximates the operation of 
the device is first derived. This circuit is then used to develop the 
computer simulation model of the PMT. Simulation results are 


presented and discussed. 
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PHOTOMULTIPLIER SIMULATION STUDIES 


I 


I . INTRODUCTION 


Simulation studies were initiated to study the effect the loading 
of the voltage-divider circuit on the photomultiplier tube response. 
Under ideal conditions the response of the PMT is linear with respect 
to light flux, and approximately exponential with respect to applied 
voltage. In practice, however, the response of such a tube depends on 
the characteristics of the circuit supplying the required 
inter-electrode voltage [1], [2]. As a result nonlinearties occur 
because of "loading" effects in the associated circuitry. An accurate 
model of the PMT is useful in choosing a suitable method of 
compensation when a very high degree of linearity is required. For 
example, the choice of the value of the capacitors on the final dynode 
stages is critical to supply peak currents for each pulse received and 
to completely recover between returns. In this report we will 
describe efforts to develop computer simulations of typical PMT tube. 
The suggested RCA 7265 base (see published data sheet) was used in the 
simulation, although our approach is quite general and can be applied 
to any type of voltage divider network used for photomultiplier tube. 

In this report we will develop mathematical models for PMT and 
PMT networks divider. Initial computer simulation are used for the 
PMT with "flat top" pulses. The light inputs are a pulse pair. A 
study is made of the effect of the first pulse on the second pulse. 
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II. EQUIVALENT CIRCUIT MODEL FOR PMT 


To simulate the PMT and its associated voltage-divider network it 
was necessary to obtain an equivalent circuit model which approximates 
the behavior of the tube under operation. To illustrate the model 
used, consider the simplest voltage divider network for a multiplier 
tube as shown in Figure 1. This model of Figure 1 is simpler than the 
RCA 7265 model, however, it is used here to explain the theory. The 
standard tube-circuited notation is used where K is the photocathode, 
Di , D 2 , . • • , are dynodes, A is the anode, Rl is the anode load 
resistance. The resistors R]^, R 2 ,... are dynodes network resistors. 
For simplicity the currents shown are "electron currents" rather than 
conventional currents which flow in opposite directions. 

From Figure 1, let cathode current be Ij. and anode current be 
I^. Assuming the currents in dynodes as shown, therefore 

10 * Ik + II (1^ 

11 ■ I 2 + Ipi (2) 

1 2 - I 3 + Id2 

In » In+1 + iDn 

lo - In+1 + Ia 

In general 

In - In+1 + iDn ’’ ^6) 


0 



In the practical model RCA 7265 shown in Figure 2 of the data 
sheet, we have 14 dynodes, i.e. n ■ 14. Also, there are capacitors 
starting at n ■ 10. 

Therefore, in this case we have 

110 * 111 + Id IdIO 

111 * Il2 + Ic2 + IdII 

112 • Il3 * Ic3 + Id12 

113 ” Il4 + Ic4 + Id13 

114 * Il5 + Ic5 + Id14 

lO • IlS + Ic5 + Ia 


Due to secondary emission [3], [4], (5], current at each dynode 
must supply a secondary emission current to the next dynode g times 
larger than the current incident to it. 

From Figure 1 the current supplied to the last dynode is: 

iDn - Ia - lA/g - Ia<1 - - Ia^-D (7) 

where r * 1/g is the ratio of the currents, i.e. the stage gain. 
Assuming, for simplicity, equal current gain ratio at each stage, we 
have in general: 

iDi * lAr^"“I^(l-r) (8) 

Ik • IaJ^" “ Ia/G (9) 

where G ■ r^ is overall gain of the multiplier. It follows from (8) 
that 

IDI • glk - Ik ' <I0> 


Let I, 


gli 


^Dl^ 


( 11 ) 
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Similarly 


'd2 ■ V 

• I,j8^(l-I/S) • Slj - 

A. 

Let - I2 

/.I 2 * ^D2'*’ ^1 (12) 

and so on for all other stages. Therefore, each successive current 
gain in each dynode stage is represented by a current source as shown 
in Figure 3. We note that in general the current gain ratio at each 
dynode stage is proportional to the voltage applied to it and it 
changes from one stage to another. 


III. CURRENT GAIN BEHAVIOR 

The secondary-emission process, which is fundamental to 
photomultipliers can be considered to consist of three physical 
processes (1) excitation of secondary electron (2) transport of these 
electrons through the solid and (3) transport though the vacuum-solid 
interface. In elementary theory of secondary emission it is assumed 
that the secondary emission yield S defined as the ratio of secondary 
to primary current is function of the primary energy [4], that is 
6 ■ j Eq) f(x) dx 

where: nCx^, Eq) is the number of secondary electrons produced at 
distance X]^ from the surface by a primary electron of energy 
Eq. f(x) is the probability that a secondary electron produced 
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at X reaches surface and is emitted in vacuum. 


Both the theoretical and experimental evidence [l], [3], [4], suggest 
that we can represent the stage gain g as power law in the 
accelerating voltage v 

.-.g - AvP (13) 

where 

A is constant depends on the material used in the tube manufacture 
p is constant generally < 1 

The constants A and p are generally determined from the tube 
manufacturer supplied curves. General discussions on the constants A 
and B are also found in [5]. For a tube with n identical stages, the 
overall tube voltage is 

V * (n+l)v (14) 

The overall gain of the tube is 

n . . p.n .n r V T” 

G . g . (AvP) . A I I 

- BV3 

where q * np or p ■ q/n, and 

q - ddog G)/d(log V) (15) 

B - A”/(n-H)”P 

For a practical tube, empirical gain-voltage curve could be used and 
hence, we can obtain q and p. Flgue 4 of RCA 7265 data sheet was 
used to determine the constants A and p. 

From Figure 4 of the data sheet, it follows;,’ 


0.8556 

(16) 

0.0448 

(17) 
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We used the above data in simulating the PMT. For the given PMT RCA 
7265 with n * 14 and from data sheet the overall gain obtained from 
the data sheet is: 

G * 4.8 X 10^ 

Therefore, the gain g per each stage is 
g - 14 Ju.8 X 1Q7 * 3.53 
Therefore, each stage has equal gain of 3.53. 

This, of course, is oversimplified as discussed above and a more 
accurate model for each stage i is 

81 * ^ 
where A and p are as given by (16), and (17) and v^ is the voltage 
across each dynode stage. We developed our algorithm based on 
Equation (18). Thus, in our results the gain per each stage 
depends on the accelerating volage of each stage. 

IV. COMPUTER SIMULATION 


The circuit shown in Figure 3 is a nonlinear circuit, which has 
nonlinear components such as the Zener diodes shown and nonlinear 
controlled sources which represents the current gain. A study of the 
linearity response of PMT and its relation to the network of resistors 
can be found in [6,7]. The component values of circuit of Figure 3 
were obtained from the data sheet of the RCA 7265 with 2400 volt power 
supply. In the simulation effort we were interested to obtain the 
variation of the anode current versus time and the variation of any 
other variable in the circuit. Each dynode gain current is 
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represented by nonlinear current source, where its value is a function 
of the voltage across it as in Equation (18). The circuit of Figure 
3, in general, is governed by nonlinear differential equations which 
are to be solved at each time step to obtain the response at various 
instants of time. Because of that, the general nonlinear simulation 
technique [8,9], may be used for simulation. In our effort, general 
equations were obtained in terms of the circuit parameters, and 
several cases were studied for different parameters. To obtain the 
cathode current, consider, for example, the two cases of the average 
photons arrival rate (it photons/sec) Ip: 
i) Ip ■ 2000 photons/pS 
ii) Ip ■ 200,000 photons/pS 

Assume quantum efficiency ■ 0.1 at 532 nm, therefore, the average 
number of photoelectrons /sec is 
Ipe ■ »llp 

For the above two cases we have 

i) Ipg ■ 200 photoelectrons/pS 
ii) Ipg » 20,000 photoelectrons /pS 
and cathode current 
Ik - Ipe * e 

where e is the electron charge ■ 1.602E - 19 coulomb 
Therefore 

1) - 3.2E-11 Amp 

ii) - 3.2E-09 Amp 
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V. 


RESULTS AND DISCUSSION 


For the following results, we assumed that the amplification 
process is noise free, that all of the current emitted from the 
cathode is collected by the first dynode, and that the seconary 
emission ratio per stage is represented by (18). We used RCA 7265 
with its voltage - divider network, and data as shown in Fig. 8 of RCA 
7265 data sheet. 

1 . Simulation of voltage variations across Rt (load resistance) with 
High Voltage Power Supply (HVPS) of - 2400 VDC with Loading Capa- 
citance Ccj ■ 500 PF 

A wavelength of 532 nm was considered. The light pulses on the 
face of the PMT of 200 to 200,000 photons per microsecond (pS) were 
used in the simulation. The light input is flat pulse pair (on-line 
pulse followed by an off-line pulse) recurring every 100 ps (10 HZ). 
Both the on- and off-line pulses were the same with a flat top, a 
width of 20 pS and a separation of 200 pS. Figure 4 shows the pulse 
applied. The capacitor C 5 * 500 PF. The pulse output across the load 
Rl is shown in Figure 5 for case of 200 photon per pS. Figures 6 and 7 
shows the pulse output across load resistance Rl for case of 2000 and 
200,000 photons per pS. In the last case the results shows 
distortion. Table 1 shows the maximum anode voltage across the load 
resistor Rl versus different average photons arrival rate Ip in the 
interval of 200 -200,000 photons per pS. A plot of this voltage 
versus Ip is shown in Figure 8 , for the case of first pulse followed 
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by second pulse. The results of Figure 8 shows that the difference 
between first and second pulse behaves as nonlinear function of 
cathode current as the light density increases. 

2 . Simulation of voltage variation across Rt with HVPS voltage - 2400 
VDC with C.^ - 1000 PF and Cg » 0.01 uF 

Next, we changed the value of the loading capacitor C 5 to 1000 PF 
and then to 0.01 pF. The other loading capacitors in the PMT network 
are changed in ratio of C 5 as in the RCA data sheet. Figure 9 shows 
the voltage across Rt (Vj^ ) for light density of 200 photon/pS and C 5 
» 1000 PF. Figure 10 shows the voltage across Rl^^Rl^ light 
density of 2000 photons/pS and C 5 * 1000 PF . Figure 11 shows the 
voltage across Rl (Vrl) for light density of 200,000 photons/pS and C 5 
* 1000 PF. Similarly, the voltage across Rj^ for different light 
density and for C 5 » 0.01 pF are shown in Figures 12, 13, and 14, 
respectively. Tables 2 and Table 3 show the maximum anode voltage 
versus different average photons arrival rate Ip in the interval of 
200 - 200,000 photons/pS for C 5 ■ 1000 PF and C 5 * 0.01 pF. A plot of 
such voltage is shown in Figure 15 and in Figure 16. 

3. Simulation of the effect of the first pulse on the second pulse. 

To determine the effect of the first pulse on the second pulse, 
the flat pulse width was increased from 20 pS to 100 pS. Figures 17, 
18, and 19 show that the voltage across Rl for different values of 
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photons/|iS and with C 5 » 500 PF. The results show that the pulse 
width variations do not affect the voltage variations across R^, where 
the maximum value is the same as obtained before. 

In all of the above simulation results it is clear that the PMT 
shows a nonlinear relation in the dynode voltages of the tube. This 
is because the secondary emission ratio depends on the voltage across 
each stage in the tube. It appears that as the number of photons/MS 
increases, and hence Ijj increases the nonlinearity increases. In 
fact, the results show that an Increase of light density more than 
4000 photons/pS causes the difference between the first and second 
pulse to behave as a nonlinear function of cathode current. All of 
the above results were simulated for flat top pulse. Further simula- 
tions for exponentially decaying pulse are needed, where a study of 
the effect of the first pulse on the second pulse should be examined. 

Additional information on PMT performance can be easily obtained 
using our developed simulation technique. For example, selecting 
optimum PMT components in the divider network and operating anode 
current for maximum linearity. Also, it is possible to use this 
technique to determine the dynode voltage change versus PMT gain. Our 
technique can be used with simple extentions in simulating PMT with 
transistor and Zener type bases or mixture of these. The above 
results are the initial phase of the simulation studies of PMT. 
Further studies as second phase of the above results are needed to 
complete the goal of the simulation. 
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VI . Recommendations for Further Study 


To fully meet the objectives of the research project, it is 
recommended that a second phase of the simulation study be performed, 
specifically to accomplish; 

1) Simulation for exponentially decaying pulse with recurring 
rates, photon arrival rates, wavelengths, turn-on delays, 
etc. Associated with the actual LaRc Electra DIAL 
observation. Also studying the PMT gain variation by 
selecting HVPS setting over the range used on the Electra 
DIAL project. 

2) Selecting optimum PMT components in the divider network for 
maximum linearity of PMT operation. 

3) Modify simulation model to include transistor base models. In 
this case simple extensions are to be added to the suggested 
simulation model. An improved base design should result from 


the simulation studies. 
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Table 1 


Maxlaim Anode Voltage for 200-200,000 Photons/jiS 


Photons Maximum Anode Voltage 


Pulse 


200 

7.68 E -03 

1 


7.68 E -03 

2 

300 

1.56 E -02 

1 


1.56 E -02 

2 

800 

3.23 E -02 

1 


3.23 E -02 

2 

1,600 

6.98 E -02 

1 


6.98 E -02 

2 

2,000 

9.09 E -02 

1 


9.09 E -02 

2 

6,000 

7.81 

1 


6.8 

2 

10,000 

2.49 

1 


3.41 

2 

12,000 

2.05 

1 


2.11 

2 

18,000 

3.73 

1 


3.64 

2 

20,000 

8.19 

1 


7.36 

2 

80,000 

1.20 

1 


11.1 

2 

100,000 

1.02 

1 


6.41 

2 

180,000 

1.41 

1 


4.33 

2 

200,000 

4.68 

1 


3.76 

2 
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Table 2 


Maxlaua Voltage Across Anode “ C5 ■ 1000 PF 


Photons Maximum Anode Voltage 


Pulse // 


200 

7.6 E -03 

1 


7.6 E -03 

2 

600 

2.3 E -02 

1 


2.3 E -02 

2 

1,200 

4.8 E -02 

1 

4.8 E -02 

2 

1,600 

6.6 E -02 

1 


6.6 E -02 

2 

2,000 

8 . 5 E -02 

1 

8.5 

8.5 E -02 

6,000 

3.5 E -01 

1 


3.5 E -01 

2 

10,000 

1.082 

1 


1.108 

2 

14,000 

3.96 

1 


4.71 

2 

18,000 

5.85 

1 


6.17 

2 

20,000 

3.146 

1 


3.543 

2 

60,000 

5.818 

1 

7.1 

2 

100,000 

14.2 

1 


30.6 

2 

140,000 

2.74 

1 

2 

3.14 

180,000 

2.50 

1 

A 


2.39 

2 

200,000 

2.621 

1 

A 


2.448 

2 
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Table 3 


Photons 

200 

600 

1,200 

1,600 

2,000 

6,000 

10,000 

14.000 

18.000 
20,000 
60,000 

100,000 

140.000 

180.000 
200,000 


Haxiaua Voltage Across Anode, C5 ” 0.01 jiF 


Maximum Anode Voltage Pulse // 

7.57 E-03 1 

7.58 E-03 2 

2.28 E-02 1 

2 . 28 E-02 2 

4.59 E-02 1 

4.61 E-02 2 

6.14 E-02 1 

6 . 18 E-02 2 

7 . 71 E-02 1 

7.77 E-02 2 

2.41 E-01 1 

2.47 E-01 2 

4.18 E-01 1 

4.40 E-01 2 

6.15 E-01 1 

6.60 E-01 2 

8.29 E-01 1 

9.16 E-01 2 

9.45 E-01 1 

1.061 2 

6.21 1 

21.04 2 

25.7 1 

28.3 2 

25.8 1 

57 2 

46 1 

3.85 2 

22 1 

57.04 2 
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UV DIAL MEASUREMENTS OF O3 PROFILES IN REGIONS OF 
SPATIALLY INHOMOGENEOUS AEROSOLS 

BY 

S. R. YOO AND R. H. LANG 


ABSTRACT 

In the DIAL measurement of the O3 concentration profile the 
corrections due to aerosol backscattering and interference in two 
wide gap wavelengths are discussed in regions of spatially 
inhomogeneous aerosols. A new Bernoulli solution for determining 
the relative distribution of aerosol backscattering is presented, 
and the boundary values for these solutions are discussed. For 
the multiple scattering effect, a factor is added in a new 
Bernoulli solution. It is shown that this factor results in a 
small correction in O3 concentration profile. 
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UV DIAL Heasxirqaents of Oj Profiles in Regions of 
Spatially Inhomogeneous Aerosols 

I. Introduction 

With the advent of Lidar systems, great advances have been made 
in recent years in the remote sensing of the atmosphere. The 
Differential ^sorption Lidar (DIAL) technique, using the Lidar 
system, has been used to obtain vertical profiles of atmospheric gas 
concentrations from ground, and airborne platforms. The major 
advantage of the DIAL technique is that the result is independent of 
many system parameters. 

This report represents a continuation of the work which has been 
done by Browell, Ismail, and Shiply in the NASA research center.^ 
They have presented the new Bernoulli solution for inverting lidar 
returns in the UV DIAL O3 measurement. In this report we reproduce 
and improve a part of the results of the new solution, and discuss 
the multiple scattering effects on the lidar returns. A brief 
summary of the new Bernoulli solution in the UV DIAL O3 measurements 
presented by Browell, et. al, is as follows. 

By using the DIAL technique the atmospheric gas concentrations 
as a function of range can be found. The measurements involve two 
types of errors, called backscattering error and extinction error. 
The backscattering error is due to interference by aerosol backsca- 
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ttering, and the extinction error is due to difference between 
Rayleigh and aerosol volume extinction coefficients at the on and off 
wavelengths, respectively. Because of the wide gap of the on and the 
off wavelengths in the UV DIAL O 3 measurements, these error terms can 
not be ignored. By using the wavelength dependency of scattering 
coefficients, the extiction error is found in terms of the scattering 
coefficients at the off wavelength, and the difference of the 
absorption cross-section at the on and the off wavelengths. Assuming 
a power relationship between the backscattering coefficient and the 
total scattering coef f icient( i.e. , power law), the backscattering 
error is expressed in terms of the backscattering ratio at the off 
wavelength, which is the ratio of the aerosol to the Rayleigh volume 
backscattering coefficients. To calculate the backscattering error, 
the backscattering ratio should be obtained. In case of aerosol 
dominated atmospheric scattering, the distribution of aerosol volume 
backscattering along the lidar path can be found by solving the 
Bernoulli equation for the lidar returns. Since the Rayleigh 
component of scattering near UV wavelength has as much importance as 
the aerosol component, the new Bernoulli equation including the 
Rayleigh and aerosol terms has been derived and solved. The quantit- 
ative analysis of two correction terms shows that the extinction 
error of Rayleigh is found to be 6.7 ppb (independent of altitude) 
and the backscatter correction ranges from - 8 ppb to 22 ppb. 
However, the extinction error of aerosol has not been implemented in 
the DIAL data reduction because of the uncertainties in the knowledge 
of aerosol extinction coefficient, and the aerosol extinction 
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correction has been remainded in the further works. 

In this paper, we try to solve the problems which have been 
suggested for further workl. One is to estimate the aerosol 
backscattering coefficient at low altitude and to obtain the gas 
concentration of interest at lower altitude from this, and the other 
is to apply the aerosol extinction correction into the DIAL O3 
measurements. In other words, the number density of gas of interest 
can be found by estimating the exact boundary value of Bernoulli 
solution at Furthermore, multiple scattering effects are 
examined to improve the inversion of O3 profile in the UV DIAL 
measurements . 

In the next chapter, the single scattering measurements are 
briefly described and in the following chapter the multiple scatter- 
ing is applied into the UV DIAL O3 measurements. By examining the 
multiple scattering effects, a small reduction of the extinction 
correction is expected. 

In the chapter IV, an algorithm for estimation of boundary value 
of the new Bernoulli solution is presented by a simple algebric 
calculation and the sensitivity of Bernoulli solution to the input 
parameters is discussed. 
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II. Single Scattering in the UV DIAL O 3 Measurements . 

Ever since the lidar system was introduced in the early 1960 's, 
algorithms for inverting lidar return signals have been discussed in 
the literatures.^'^' Until recently, however, the inversion of lidar 
returns has been plagued by instability and inaccuracies. The 
instability of previous solutions is found to be caused by mathema- 
tical and not, as some have suggested, physical reasons. Inaccura- 
cies can occur by making unnecessary assumptions to further simplify 
the equation. The solution proposed by Klett® in 1981 was a step to 
improve this situation. 

A. Lidar Equation . 

The lidar equation is the governing relation which permits the 
calculation of the extinction coefficient and estimates of concentra- 
tion under conditions of single scattering from a given lidar 
return. The lidar equation can be derived from the radar equation, 
the equation of radiative transfer, or the equations related on the 
laser radiation. It is convenient to begin with the equation of 
radiative transfer since it reveals the physical assumptions made in 
the lidar equation: 


dIf(R,s) 

dR 


= - *If + 


Ct(R) 


4n 


P(s,s' ) *If (R,s) ‘dfl’ + e(R,s) 


4n 


= - a(R) • [If (R,s) - Jf(R,s)] 


( 1 ) 


where If(R,s) is the specific intensity of radiation at frequency f 
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in the s direction, a(R) is the total scattering and extinction 
coefficient, and Jf{R,s) is the source term. The source term 
contains contribution to scattering from all s’ direction and 
emission in the s direction, i.e., e(R,s). P(s,s’) is the scattering 
phase function. In the radiative transfer equation, an important 
assumption has already been made: that the distance between the 
scatterers is larger than the wavelength, and the scattered radiation 
is incoherent and interference effects are not considered. 

For the next step in deriving the lidar equation, the assump- 
tion is made that If(R,s) >> Jf(R,s). This assumes that there are no 
significant amount of emission or scattering from the medium in the 
direction of propagation. Note that this contains the critical 
assumption that only single scattering events occur. In making this 
assumption we obtain 


or 


1 dlf(R) 
If dR 


- 


If(R) 


I© exp{- 


af{r' )dr’ } 


( 2 ) 


(3) 


with Iq being the initial intensity of the beam before any extinc- 
tion. We note that the equation (2) is a function of range only and 
equation (3) is the well-known Beer’s law. 

For the lidar system in which the return power is being received 
in the direction opposite to the beam direction, the return power 
will be proportional to the initial transmitted power Pq, geometric 
attenuation A/r 2, where A is the receiver area, and the folded laser 
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pulse length given by cx/2, where c is speed of light and x is the 
time duration of the laser pulse. Furthermore, the return power will 
depend on the volume backscattering coefficient, P|^(R), at range R, 
and finally on a system function F(R) which includes parameters such 
as system sensitivity and geometric crossover of the receiver and 
laser beam. Putting all these factors together with equation (3), we 
obtain the lidar equation: 


CT A 

P(R) = Pq — F(R)•^TI(R) — exp{ -2 


fR 


r2 


a{r) dr} 


(4) 


where the extinction term in Beer's law has been squared since the 
return power must traverse the same medium twice. 

As we discussed above, the assumptions in the lidar equation (4) 
can be summarized as followings: 

1. incoherent scatterers (Equation of radiative transfer), 

2. single scattering. If >> Jf, 

3. only backscattering is considered, 3y^(R), 

4. Monochromatic ligt source. 

5. Beam overlay is complete, i.e., no blind spot. 
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B. Differential Absorption Lidar (DIAL) Technique. 

Differential Absorption Lidar is a single ended active remote 
sensing technique for monitoring concentration of various atmospheric 
gases. In the DIAL technique two laser pulses of slightly different 
wavelength, where one is chosen to coincide with a strong absorption 
feature of the specific constituent of interest and the other lies in 
the wing of this feature, are used and the attenuation of the two 
pulses is compared. The wavelengths of two pulses will be called the 
on and off wavelengths in this paper. 

For a monostatic single-frequency pulsed lidar, the single 
scattering lidar equation is obtained in equation ( 4 ) . The ratio of 
the return power signals at the on and off wavelengths, i.e. ,P(X q^/R) 
and P(XQff,R), respectively, is expressed in this form: 

R(Xon^^) Xq|.j,R) ^ 

= exp[-2 {a(Xon^r) " a(Xoff }dr] (5) 

R(Xoff,R) Pxx^^off'^^ 

where the subscripts of on and off refer to the on and off wave- 
lengths. In equation (5) we assume that the all system parameters 
at both wavelengths are equal: 

F(Xon^R) = F(Xoff,R). (6) 

The total volume extinction coefficient, a(X,r), is composed of 
scattering terms and absorption terms of atmospheric components 
(e.g., gases and aerosols). Hence: 

a(X,R) = ag(X,R) +aa(^,R). (7) 

In the DIAL technique, the absorption of the lidar power is 
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mainly caused by the gases of interest. The total volume extinction 
coefficient of absorption, then, is the sum of the extinction coeffi- 
cient of absorption due to each gases; 

cta(A/R) = ? aa(A,R) =N(R)-Oa(A) (8) 

where N(R) represents the number density of the gas of interest which 
causes the absorption, and o^ is its absorption cross section. 
Substituting equations (7) and (8) into (5), and taking log and 
integral yield: 


1 P(Aoff»R) 1 pR^Acff'^^ 

N(r)*Aoa*dr = - In - - In 

o 2 P( Aqjj,R) 2 Aon^^^ 


^s^^on'^^ ~ ^ ^of f ' ^ ^ 


where we have introduced the differential absorption cross section, 
Ao^, that is the difference of the absorption cross section at the on 
and off wavelengths: 


^^a = «^a(^on) " <^a<^off> (10) 

If we define the average number density of gas and assume that 
the volume scattering coefficient ag{r) is constant over the range 
increment (interval), AR = R 2 - Ri, i.e., 

1 PRt 

N » N(r)-dr (11) 

^1,^2 


and 
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ctg( A) 

Ri,R 2 


1 [R 2 

as(A,r) 

AR jRi 


AR 


(12) 


where Tq represents the differential extinction spatial depth for the 
atmosphere . 

From the assumption of (11) and (12), the range information of 
the gas of interest can be obtained in the interval, R^ and R 2 : 

1 P(Aon'^l^ P(^off»^2^ 

N = In {M} 

Ri,R 2 2AR«Aog^ ^^^on'^2^ P(AQff,R]_) 


1 ^Ti^^on'^l^ I^Tt^ Aoff /^i) 

In 

2 A R'zia-j ^n^^on'^2^ ^n^^off'^2^ 


1 

TT ^®^s^^on^ " ®'S^AQff)} 

2AOa Ri,R2 


{B} 


{E} 


(13) 


For the O 3 number density in the DIAL measurements, we rewrite 
equation (13) in this form; 


N( 03 ) = {MO 3 } + {BO 3 } + {EO 3 } (14) 
Ri,R2 

The terms {BO 3 } and {EO 3 } in equation (13) are the corrections 
to the DIAL measurements due to the wavelength dependence of bac)c- 
scattering and volume scattering of atmospheric gases. Under condi- 
tions of spatially homogeneous atmosphere the terra {BO 3 } can be 
neglected, i.e., since “ 3 ti(R 2^* Also, when the on and off 
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wavelengths are close together, the term {EO3} generally becomes 
negligible. Then only term {MO3} can be used to evaluate the range 
concentration in the DIAL measurements. 

In the UV DIAL measurements for the density of O3 , neither of 
these conditions are present under normal atmospheric circumstances 
because ozone is usually distributed inhomogeneous ly in the normal 
atmosphere and the wide gap between wavelengths for the on and off 
frequencies are used in the O3 measurement (286 nm and 298.3 nm) . 
The error terms, i.e., {BO3} and {EO3 ), will be regarded to as the 
backseat ter error and the extinction error, respectively. These two 
terms for the O3 measurement will be discussed in the next section. 
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C. Two Interference Terms in the O3 Measurements . 

In most common scattering situations, separation of the gaseous 
(Rayleigh scattering) and aerosol (Mie scattering) components from 
the lidar measurements is needed. The volume backscattering coeffi- 
cient and the volume extinction coefficient due to scattering can be 
expressed in the svim of their coefficients of the gaseous and aerosol 
components : 


X,R) ~ Pn,R^ X/R) t 3xt,M^ X»R) 


where subscripts R and M represent the gaseous (Rayleigh scattering) 
and aerosol (Mie scattering) components, respectively. 

The volume scattering coefficients for gaseous and molecular 
components are dependent on the wavelength. These coefficients vary 
inversely with wavelength as powers, 4 and 6, respectively (Collis and 
Russell^) : 



-4 

R(X) 

* <=^s,R*^ 


-6 

m(>^) 



(16) 


where and 0 , 5 , M average scattering coefficients indepen- 

dent of wavelength in range R^ and R 2 . 

Using Taylor expansion, i.e., f(x) * f(Xo) + f'(Xo)(x - Xq)/ we 
can write the volume scattering coefficient at the on wavelength by: 


®’S^^on^ “ ®s^^off^ ®*s^^off^*^ ^on ” ^off^ 


(17) 
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or defining that AA = Aqjj - A off / we find; 

Q-s^^on^ “ ®s^Aoff^ ® ®'s^Aoff^' 

= ■'■ ( 4 -aR(Aoff) + 6*OM(Aoff)}. ( 18 ) 

^off 

Therefore, the extinction error {EO3} in interval range and 
R2 in equation ( 14 ) can be approximated in this form; 

1 

{EO3) { 4 *ag^R(Aoff) + 6*^5 M^^off^ ( 19 ) 

Ri,R 2 Og Aoff 

The scattering coefficients at the off wavelength do not cause 
absorption. These coefficients can be obtained from the relation; 

ctg ^R(A,R) s* Nr( ® ^ *^S , r(^ ^ * 

( 20 ) 

where N|^(R) and Nj,j(R) represent the number density of gases (Rayleigh) 
and aerosols, and Og^R(A) and cjg^j5(A ) are their cross section. 

Typically the Rayleigh number density can be ta)cen from the 
U.S. standard atmospheric number density. 

The backscattering correction {BO3} can be evaluated by intro- 
ducing the backscatter ratio that is the ratio ' of the backscattering 
coefficients for gaseous and molecular components; 


b^(R) 




Pn,R(l/R) 


( 21 ) 
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In equation (21), the relation of volume backscattering 

coefficient and volume extinction coefficient due to scattering is 
needed to further simplify the backscatter ratio. The traditional 
way to account for any dependence between the and Cg is to assume 
a power relationship of this form: 

k 

,R) = BQ*ag(A,R) (22) 

where k depends on the lidar wavelength and various properties of the 
obscuring aerosol, and Bq is constant. Klett^ has indicated that the 
value of the exponent k is generally in the range 0.67 S k < l.o. 
Although k is no longer unity, the ratio of backscattering to 
extinction coefficients, Bq, can be found by assuming that k = 1 in 
absence of absorption. We will consider later the case of k = 1. 
The backscattering/extinction ratios for Rayleigh to aerosols 
components in general have been expressed in the form: 


0ti,r{^ /R) 


3 

- — ®s,R<^ 




(23) 


where = Pj^/4n is the normal aerosol backscattering phase function 
is taken to be a constant for the isotropiq* scattering. 

Noting the wavelength dependence in the extinction coefficients 
,i.e.. Eg. (16), the backscatter ratio can be expressed in the form:^ 


b,(R) 


-6 

-4 

3/8rt*ag^I^(R) 


4-6 

bo(R) * A 


(24) 
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The backscatter ratio at the on wavelength in terms of that at 
the off wavelength becomes: 

Aon 


off 


= boff (R) • (1 - 


Al. 


4-6 


off 


Knowing that A A = Agff ” A qj^ is small, we have; 
bon(R) * boff(R).{l - (4-6)— A - ^ -) 


off 


or 


AA 


bon(R) - boff(R) * (4-6)—- boff(R) 


off 


(25) 


(26) 


(27) 


To obtain the {BO 3 } in terms of the backscatter ratio, we 
consider that: 


{ } 


^on i " ^off(^) 

= ( ) • {1 + > (28) 

^off 1 + ^3off<R) 


where the wavelength dependence in the backscattering coefficient of 
Rayleigh component at the on and off wavelength has been assumed. 

By substituting (27) into (28) and taking logarithm, we have: 


ln(l - 4- 


AA 


•) + ln(l 


^off 


(4 - 6 ) • AA boff (R) 

A off 1 + boff(R) 
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{ 4 - 6) • AA ^off 

s - 4—^ (29) 

Aoff Xoff 1 + boff(R) 

where the approximation in equation (29) has been made, i.e., 
log(l+x) ® X, as X << 1. 

Finally, the baclcscatter correction {BO3) in interval range 
and R2 in equation (14) is obtained in the form; 

(4 - 6) AA boff(R2) ^off(^^l> , 

{BO3} * { - } (30) 

2* R* Oa Xoff 1 + boff(R2> 1 + ^^off^^l^ 


To calculate the {BO3}, we need to evaluate the backscatter 
ratio at the off wavelength. In other words, the aerosol 
backscatter ing profile along the lidar return path must be obtained. 
From equations (21) and (23), we can rewrite the backscatter ratio in 


terms of the function g,(R) by: 


b,(R) = 


g,(R) 


(31) 


where the function g (R) is represented by: 


3/8n*as,R,( A,R) 

g,(R) = 

^ Pn(R) 


(32) 


Since the Rayleigh scattering coefficient can be obtained in the 
relation in equation (20), g^(R) is a known function. As mentioned 
above, the aerosol distribution profile along the lidar path needs to 
be found. It can be obtained by solving the Bernoulli equation which 
contains the Rayleigh and the aerosol scattering components. 
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D. Inversion Methods and Hew Bernoulli Solution for the UV DIAL 
O 3 Bteasurements . 

From the single scattering lidar equation, the inversion of 
extinction coefficients has been discussed in the slope and ratio 
method, Bernoulli method ( forward or backward inversion methods), AGIL 
method^, etc. In these methods, the basic relation of equation (22) 
has been assumed. Noting that the following definition is made: 

S(R) = ln[R2.p(R)] (33) 

we can derive the differential equation for the extinction coeffi- 
cient: 

ds k da 

— 2a (34) 

dR a dR 

In deriving equation (34), the system parameter F(R) and Bq are 
assumed to be constants. 

The following sections briefly review the major inversion 
methods and the new Bernoulli solution for the UV DIAL O 3 measure- 
ments. 

1. Inversion Method of the Lidar Equation . 

The brief review of inversion of Lidar return signals is 
summarized in the followings: 

Slope and Ratio method 

If one assumes that over some distance the aerosol is uniform, 
i . e . , da/ dR = 0 . , we can have : 
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1 dS 

a (35) 

2 dR 

so that the slope of a plot of S versus R will give us the value of a 
over the homogeneous region. This can be however a very poor 
approximation if there is noise or if the medium is inhomogeneous. 
The ratio method is just a variant of the slope method when it is 
applied segment by segment. Note that the equation would result from 
the assumption that k =0, i.e., the backscattering is a constant 

independent of the aerosol concentration. The slope method has been 
used most often in the past since it is simple. 


Bernoulli method . 

The equation (34) can be solved exactly because it is a standard 
differential equation, ncunely the Bernoulli equation or a homogeneous 
Ricatti equation. The solution, which is called the forward integral 
method, is; 


a(R) 


exp[ (S - So)/k] 


-1 


{OIq 



'R 

exp[(S - SQ)/k]dr> 


(36) 


where the subscript o refers to the reference distance Rq which is 
the closest to the lidar system. This solution is however unstable 
since two relatively large terms are subtracted in the denominator to 
obtain a small difference (Klett, 1981)®. Thus, if there is an error 
in Oj^, the solution quickly becomes unstable, and an unrealistic 


18 


value of a{R) are computed at modestly small R. In order to obtain 
the stable solution of a(R), the idea of making a guess or measure- 
ment of a(R) at the largest distance from the lidar system, i.e., Oj^, 
was mentioned in the reference. However, after some unsuccessful 

attempts, Klett obtained the following solution which is called the 
backward integral method:® 


a(R) 


exp[(S - S„,)/k] 


-1 2 
tOm + — 
k 


exp[ (S - Sjn)/k]dr> 
R 


(37) 


where the subscript m refers to the reference distance Rj^ which is 
the furthest from the lidar system. Now that the two expressions 
were added in the denominator, the solution becomes more stable. 
Furthermore, the error in a„, becomes less important as R decreases 
and the integral in (37) increases. Thus the solution will tend to 
converge to the correct a(R) as R gets smaller. However, the 
boundary value 0 ,^ must still be determined. 


2. New Bernoulli Solution in the UV DIAL Ot Meastirements . 

In order to construct an Bernoulli eguation in the DIAL techni- 
que we define the normal lidar return in the form; 


P(X,R)-R^ 

Ux(R) = r 

^ P(\,Ro)-R^ 

where the subscript o refers to the reference 
the closest to the lidar system. 

Recalling the lidar equation (4) and 


(38) 

distance R© which is 
the relations of the 
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^ fR 

U>(R) = *exp{ -2 ta-(X,r) + N(r) *aa( X) ]dr} (39) 

^ 3n(X,Ro) jRo 


The ratio of the backscattering coefficients in range R and Rq 
can be written in terms of the backscatter ratio defined in (21): 


Pn( X .R) 0 ti,r( X ,R) 1 + b, (R) 

= — : { i } (40) 

X»Rq^ 3ti,R^^»Ro^ ^ bj^(RQ) 


By defining that the scatter ratio is: 


»Ti( ^ /R) 

4>,(R) = = 1 + b,(R) (41) 

Pti,r( ^ /R) 


and from equations (15) and (20), we find: 


Nr(R) (R) fR 

u (R) = •exp{-2 [Cg R(X,r) + Cg N(r) •Ca(X ) ]dr) 

Nr(Ro) ♦x^Ro) JRo ' 

^ (42) 

where Nj^(R) and N(R) refer to the Rayleigh gas numlaer density and the 
gas number density absorbing the lidar power, respectively. 

Since the Rayleigh component of scattering near UV wavelength 
has as much importance as the aerosol component, the extinction 
coefficient due to Rayleigh component is specified in equation (42). 
Let ^(R) represent the normalized lidar return corrected for 
Rayleigh scattering and eOasorption such that: 


Np(Rq) fR 

uJ*(R) = u,(R) *exp{ +2 [cis,Rt X/J^) + N(r) *ag(x ) ]dr} (43) 

Njj(R) JRq 
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Substituting (43) into (42), we get; 


^ <*• (R) 

U*(R) = exp{ -2 

^ ♦.(Ro) 


,r) -dr} 


( 44 ) 


where the absorption term in (42) and (43) can be neglected at the 
off wavelength. The extinction coefficient due to the aerosol 
scattering, written in the expression; 


ag,M<A.R) * g^(R)*b^(R) = g^(R)-[*^(R) - 1] (45) 

where gj^(r) and b^(r) have been defined in the equations (21) and 
(32). Then taking log, differentiating, and rearranging the equation 
(44), we have the form of a differential equation; 


d* d 

— = *(R){ — In u*j(R) - 2*g(R)) + 2'g(R)**2(R) ( 4 g) 

dR dR 


This is the inhomogeneous linear equation of Bernoulli equation; 

y’ = a(x)*y + b(r)-y2 (47) 

with the solution; 

Cl 1 fx 

y = b(t)*I(t)*dt ( 48 ) 

Kx) l(x) J 


where 


rx 


Kx) = exp( 


a(t) *dt) 


(49) 


Following the steps above, we have found the solution of the 
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Following the steps above, we have found the solution of the 
equation (46); 


U*(R) *exp[ 

X 


♦, 11 .) . 


Jr« 


r)dr] 


-1 

♦o - 2' 


dr [g(r)«u*(r) ‘exp{ 
Ro ^ 


‘r 

-2 g(i 

jRo 


(50) 


r’ )-dr'}] 


where u*(Rq) = !•/ and is the boundary value at range Rq which is 
the nearest to the lidar system. 

Using the idea of the baclcward integration method, we can obtain 
the stable solution of the equation (46) in the expression: 


* 

u (R)*exp[-2 g(r)dr) 

JR 

♦ (R) * 

^ fRm * fRm 

+ 2«j dr[g(r) (r) •exp{-2j g(r')*dr’}] 


(51) 


where is the boundary value at range R(„ which is the furthest to 
the lidar system. 
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III. Multiple scattering In the UV DIAL 0-> Measurements 

In the previous chapter, the single scattering effect was 
considered in the inversion of the lidar equation. In this chapter, 
we discuss the multiple scattering effects on the lidar equation and 
the UV DIAL O3 measurements. 

A. Multiple Scattering corrections in the Lidar Equation and the UV 
DIAL O 3 Measurements . 

As mentioned earlier, the multiple scattering effects have been 
considered in the inhomogenious or dense media such as cloud, fog and 
rain falls. When the cloud is dense enough, significant multiple 
scattering will invalidate the lidar equation. Thus, either a new 
lidar equation incorporating multiple scattering should be used or 
the data adjusted so that only the single scattering component of the 
signal is considered. 

In the earlier papers, a simple expression of lidar 
equation with multiple scattering factor has been obtained in the 
form; 


CT A 

P(R) = Pq F(R)-Pti(R)— r*exp{ 

2 R^ 


fR 


-2 


[a(r)*(l 


r(r))]dr} 


( 52 ) 


where we note that the multiple scattering effect, rir) is taken into 
account inside the exponential term. 

By examination of the available theories and experiments, 
equation ( 52 ) would be considered when the optical depth is greater 
than 0.1 for lidar systems approximately 1 Ion or more from the media 
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and perhaps for optical depths greater than 1.0 when the lidar system 
distance is close (about 50 - 100 m) . If the contribution of 
multiple scattering is small, the multiple scattering factor can be 
taken to be small and a constant. 

In the DIAL technique, we can break the extinction coefficients 
into two terms due to scattering and absorption. Since the multiple 
scattering does not dominate the absorption, equation (52) can be 
written in the expression; 


CT A 

P(R) = Pq F(R) •- ^*exp{-2 

2 r2 


[as(r)-(l - r(r)) + aa(r)]dr> (53) 


Assuming that the multiple scattering factor is small and a 
constant, and from the equation (8), we can write the equation (13) 
in the form : 


1 P(Aon»^^l) P(Aoff'^2^ 

N = In —————— 

Ri,R 2 2'AR-AOa ^(^oti'^2^ 


1 ^Tt^^on'^l^ ^u^Aoff'^2^ 

- In 

2*AR*AOg^ 0Tt^^off,Rl) I^TX^^on'^2^ 

1 

“ " r’)*{as(Aon) “ cis(Aoff)>' 

^ Ri,R2 


( 54 ) 
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B. Corrected Bernoulli Solution In the UV DIAL Ot Measurements , 

As defined earlier, the normalized lidar return, taking the 
multiple scattering efffects into account, becomes: 


— ^Ti ^ X » R ) 

Ux(R) “ ‘exp{-2 

^Tl ^ X / ^ 


fR 


[as(X,r)-(l - r) + N(r)-aa(x)]dr> (55) 


From equation (42), and by neglecting the small contribution of 
multiple scattering in the Rayleigth component, equation (42) becomes 


Ux(R) = 


Nr(R) *(R) 


Nr(Rq) *(Ro) 


-•exp{-2 


fR 

- r) + as,R(X,r) 


+ <56) 

where Np(R) and M(R) refer to the Rayleigh gas number density and the 
gas number density edasorbing the lidar power as defined earlier, 
respectively. (R) is the scattering ratio as in equation (41). 

A 

The normalized lidar return corrected for Rayleigh scattering 
and absorption is defined in equation (43): 


— * 


Nr(R) 

uHR) * Uii{R)‘ *exp{ +2 

Nr(Ro) 


fR 

<®^s»R<X'J^) + N(R)-aa(x )]<ir} (57) 

Ro 


Substituting (57) into (56), we find: 


u*(R) 


4 (R) fR 

exp{-2 as,M<X/r)-(l - D-dr) 

♦;^(Ro) JRo 


(58) 
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Following the steps in equations from (45) to (50), we found the 
stable solution of the corrected Bernoulli equation in the expre- 
ssion: 


U*(R) •expI-2 


*^(R) = 


^ (r)'dr] 
R X 


2 




(59) 


dr[g^ (r) *W^(r) *exp{-2 


R_ 

g (r' ) -dr' }] 
r ^ 


where g(r) is defined in equation (32) with )c - 1. 


g,(R) 


3/8Ti-as,R(X ,R) 



( 60 ) 


where k can not be taken to be k s 1. 



26 


IV. Analysis of DIAL Data for the O -^ measurement of in the OV 

Wavelength 

DIAL data for the O 3 measurement in the UV was obtained from the 
on-board LIDAR system on 7 Aug. 1980 (Figs. A, B). In order to invert 
the Lidar return signal for the O 3 profile the new Bernoulli equation 
was derived and solved in the preview chapter. The solution of the 
new Bernuolli equation, i.e., 1 + b^(R), is given in Eqs.(50) 

and (51). For the reasons of numerical stability, backward integra- 
tion of the Bernuolli solution is preferred to integration in the 
forward direction. However, the boudary condition at range Rmfi.e*/ 
is unkown. Recently several of papers which describ the 
algorithm for estimating the boundary condition of the Bernuolli 
solution has been published by Ferguson and Stephens,^ and Mulders^^ 
The algorithms of Ferguson and Stephens solved this problem by 
iteration, and the algorithm of Mulders solved the problem by simple 
algebra. Now we find the boundary value at range Rmr for the new 

Bernoulli solution by simple calculation in the following section A. 
In addition, some comments are made concerning the new Bernoulli 
solution at the low altitude. In the section B, we discuss the 
sensitivity of parameters such as k, and in the section B. 
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A. Bstlmatlon of ta. Boundary Value. 


The new Bernoulli solution requiring backward integration from 
some range to R has been found in the Eg. (51). 


(R) = 


i*(R) «exp[-2 •dr] 


+ 2 


'Em 

[g(r) •u*(r)*exp{-2 g(r’ ) *dr ' } ]dr 
R ^ Jr 


Note that the normalized lidar return and the corrected lidar 
return at the off wavelength has been defined in the forms. 


Uoff(R) = 


P(R) -r2 


(61) 


Neglecting the absorption term at the off wavelength, the 
normalized lidar return for Rayleigh scattering extinction becomes. 


and 


NR<Rin) 

U (R) = Uoff(R) — — — *exp{ +2 

Nr(R) 


3/8n*as r(R) 
goff(R) « —If 


■Rm 

ag,n(r)dr) 

R 


(62) 


The Rayleigh scattering extinction coefficient at wavelength,^ , 
and altitude Z, is defined by 


®s»R^^» 2) = P{((Z)‘ Oj^(x) 


(63) 
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where the total Rayleigh scattering cross section Or(X), and the 
particle number density of air Pj^(Z) are assumed to be( Measures, 
1984) 


OrCX) 


8n 

3 


6 + 36 


On X { 


P 


6 - 76, 


[cm^sr^ ] 


(64) 


550 4 - 28 

o„(x) * 5.45 X { } X 10 [cm^sr”^ ] 

X(nm) 


Pr(Z) = N* exp{ -Z/9.5 } 


[molecules cm“^] 


Assume that the depolarization ratio 6p for air is taken to be 
0.042 and the particle number density of sea level, N, is taken to be 
2.55 X 10 [molecules cm ] from the U.S. standard atmosphere. Then, 
the scattering extiction coefficient at off wavelength 298.3 nm and 
altitude Z, can be found in the expression, 

-1 -1 

®s»r( 2) = 0.1444567 exp{ -Z/9.5 } [ km sr ] (65) 

Let Z^ distance of the Lidar system from the ground, and R 
distance from the Lidar system to the actual altitude to be measured. 
Then the actual altitude beccmies Z * Z^ - R or R = Z^ - Z, Inspect- 
ing and applying the above equations, we find the expression. 


z-z, 


m. 


u(Z«- Z)*exp{ )*exp[Ci* {ex; 

9.5 



-Z 


m 


- exp—}] 
9.5 


-1 

♦tn + C2 




2 ra - “2 

dz' -utZa-z’ ) •exp[Ci{exp— 
Z 9.5 



♦off(Z) * 


( 66 ) 
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where the constants Cl and C2 , are obtained by 

3/8n 

Cl = 2.746773 X { 1 } 


3/8n 

Ct = 0.2889134 X X exp{-2ra/9.5> 

Ptt 


( 67 ) 


The solution of the new Bernoulli equation was defined in the 


form. 


&Tl^^Off 

♦off(Z) = = 1 + boff(Z) 

Pn,R^^off 


( 68 ) 


and the backseat ter ratio at off wavelength was defined by 


boff(Z) = 


^n,M^^off ®s,M^^off'^^ 
Pn,R^^off » 2 ) 9 off^ 2 ) 


The backscatter ratio at the highest altitude , Zq , ( i . e . , over 
2000 meter) is assumed to be zero and then we can express the 
Eg. (51) in the form; 


♦off^^o) = 


U(Z.- Zo)*exp{ }-exp[Ci*{exp — ^ exp }] 

9.5 9-. 5 9.5 

^m — 1 

♦m ^2* dz' •u(Za-z' ) *exp[Ci{exp exp }] 

. Zq 9.5 9.5 


(69) 


where 


for Zq > 2000 meter 
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Hence, the integral constant or the bountary constant, 
altitude Z„, can be found in the expression; 


-1 

= u(Za- Zo)*exp{- 


] 

9.5 


fZ, 


- C2- 


m. 


dz' •u(Za-2' ) *exp[Ci{exp 


-Zo 


-Zm 

9.5 


9.5 

-z’ 

- exp 

-Zm 

9.5 

9.5 


}] 


(71) 


The Eg. (71) shows that one may obtain the correct using scalar 
iteration for the small deviation of the backscattering ratio at the 
highest altitude Zq. However, the results by the slope method 
indicate that ^off^^o) is near 1.0 and it makes sense physically at 
high altitude over 2000 meter. This boundary value can now be 
used in the new Bernoulli solution requiring forward integral, 
i.e. ,Eq. ( 50 ) , for the lower altitude(i.e. , Z < Zjjj) . This is accomp- 
lished by replacing z with Z„j and performing the required integra- 


tion. 
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B. Sensitivity of the New Bernoulli Solution to Input Parameters. 

In this section, we examine the sensitivity of the Bernoulli 
solution to input parameters from the data set, (Fig., A,B), provided 
by NASA, Langley Research Center. The $(z) is computed by using Eq. 
(69) where the boundary value, 4^, is determined in the equation (71). 
In the equation (71), the boundary value is found in the range from 
1.5 to 2.8 depending on the lowest altitude. Since the data provided 
by NASA is not accurate in the lower altitude range, three values of 
have been chosen in order to complete the calculation. They are 
= 2.0, 2.4, 2.8. Figure C shows the profiles which have been obtained 
for each case of the above three choices of Figures D and E show 

that the backscattering correction varies from -20 ppm at 1000 meter 
to 40 ppm at 2000 meter. The extinction error also varies from 6 to 
20 ppm. This result includes a multiple scattering factor of 0.2. 
The corrected 0^ profile is shown in the figure F. Comparing the 
uncorrected 0^ profile, it reduces the sharpness at 1500 meter. 


I 
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FIG. A THE RETURN POWER SIGNAL AT THE OFF WAVELENGTH 





RATIO OF THE RETURN POWER 
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FIG B RATIO OF THE RETURN POWER AT ON AND OFF WAVELENGTH 




FIG. C BERNOULLI SOLUTION AT THE OFF WAVELENGTH (298.3) 




PPB CONCENTRATION 


FIG D BACKSCATTERING CORRECTION - PPB 




FiG E EXTINCTION ERROR PROFILE - PPB 




PPB CONCENTRATION 

FIG F UNCORRECTED AND CORRECTED OZONE PROFILE 



